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ABSTRACT 


The  Birnbaum-Saunders  fatigue  life  distribution  with  shape  para 
meter  o  and  scale  parameter  £  is  considered.  The  scale  parameter  Bjis 
also  the  median  lifetime,  and  assuming  that(B  is  known,  Bayes  estimators 
of  the  reliability  function  are  obtained  for  a  family  of  proper  conjugate 
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priors  as  well  as  for  Jeffreys’  vague  prior  for (oj4  When  both  (a)  andJS; 
are  unknown,  a  modified  Bayes  eat imator  ^(t^  of  the^ reliability 
proposed  using  a  moment  estimator  of^jl.  In  addition  to  being  comp¬ 
utationally  simpler  than  the  maximum  likelihood  estimator  (mle)  of 
reliability,  Monte  Carlo  simulations  for  small  samples  show  that  Rj[(tV 
is  better  than  the  method  of  moments  estimator  for  all (5  and  as  good 
as  the  mle  for  small  o  in  the  sense  of  mean  squared  errors.  _ 


Key  Words :  Patigue  life  data;  Bayes  estimation;  Method  of  moments; 
Maxima  likelihood* 


1 .  INTRODUCTION 


With  the  relatively  small  amounts  of  fatigue  life  data  that  can  be 
obtained  for  a  specific  type  of  structure,  many  of  the  two  parameter 
lifetime  models  such  as  the  gamma,  lognormal,  or  Weibull  distribution  can 
be  fitted  reasonably  well  with  respect  to  mean  end  variance.  However, 
there  is  a  wide  discrepancy  among  the  models  for  predicting  the  higher 
percentiles  of  the  fatigue  lifetime  variable.  Blrnbaum  and  Saunders 
(1969a,  1969b)  proposed  a  two  parameter  life  distribution  for  material 
failure  due  to  fatigue  crack  extension  under  cyclic  loading. 

The  Blrnbaum  -  Saunders  fatigue  life  distribution  has  probability 
density  function 

f (t;  a,B)  -  (2^57  o  6  t2)“l  (t2  -  62>/C(t/S>*  -  (B/t)**] 

x  exp  [-(t/B  +  B/t  -  2)/2o2],  t  >  0  ,  (1.1) 

with  o,B  >  0  and  reliability  function 

R(t;o,B)  -  «  [-((t/B)**  -  (B/t^/ol,  t  >  0  ,  (1.2) 

where  ♦(• )  denotes  the  standard  normal  cumulative  distribution 

function.  The  parameters  o  and  B  can  be  interpreted  as  shape  and 

scale  parameters,  respectively.  Also,  B  is  the  median  life  (Saunders, 

2 

1974)  and  the  mean  and  variance  are,  respectively,  v  ■  6(1  +  a  /2) 

2  2  2 

and  o  ■  (aB)  (1  +  5a  / 4)  .  Other  properties  of  the  distribution  are 
described  by  Blrnbaum  and  Saunders  (1969a,  1969b)  or  Mann,  Schafer, 
and  Singpurwalla  (1974).  Saunders  (1974)  studied  a  £ -normal  family  of 


2 


random  variables  closed  under  reciprocation,  of  which  the  density  (1.1) 
is  a  member.  He  discussed  moment  estimation  as  well  as  other  estimators 
of  6  ,  since  the  maximum  likelihood  estimator  (mle)  of  0  cannot 
be  found  in  closed  form. 

In  this  paper  we  consider  estimation  of  the  reliability  function 
(1.2)  by  Bayesian  and  "modified  Bayesian"  techniques.  Bayes  estimators 
for  a  proper  conjugate  family  of  priors  on  o  and  vague  priors  on  a  will 
be  given  in  Section  3,  assuming  0  is  known.  This  is  reasonable  since  0 
is.  the  median  or  "characteristic  life,"  and  information  concerning  0 
may  be  known.  In  Section  4  a  modified  Bayes  estimator  of  (1.2)  is  proposed 
and  compared  with  the  mle  and  method  of  moments  estimator  (mme)  of 
(1.2)  for  small  samples.  Monte  Carlo  simulation  results  Indicate  that 
the  modified  Bayes  estimator  Is  better  than  the  mme  but  not  uniformly 
as  good  as  the  mle  .  However,  the  mle  must  be  computed  by  Newton's 
or  other  numerical  iteration  procedures,  whereas  the  modified  Bayes  and 
mme  are  easily  computed  in  closed  form.  An  example  is  given  in  Section  5. 

2.  ESTIMATION  OF  a  AND  0 

Point  estimation  of  o  and  0  has  been  investigated  by  Bimbaum 
and  Saunders  (1969b)  and,  for  the  more  general  £ -normal  family,  by 
Saunders  (1974).  Their  main  results  will  be  stated  briefly  in  this 
section  for  completeness. 

Let  X,,  ...,  X  denote  a  random  sample  from  the  pdf  (1.1).  Birn- 
i  n 

A 

baum  and  Saunders  (1969b)  showed  that  the  maximum  likelihood  estimate  0 

2 

of  0  la  the  unique  positive  solution  of  g(x)  ■  x  -  x[2r  +  K(x)3  -t 
r[s  +  K(x))  ■  0  ,  where 


3 


s 


» 


and 


K(x) 


(x  +  X^-1  ]_1 


are  the  harmonic  mean,  arithmetic  mean,  and  harmonic  mean  function, 

A 

respectively.  In  addition,  r  <  B  <  s  .  The  mle  of  a  is  then 

A  A  A  t 

a  *  (s/8  +  8/r  -  2)^  and  the  mle  of  reliability  (1.2)  is  given  by 

R(t  ;  a,  3)  »  ♦E-«t/6)*s  -  (B/tjS/a]  . 

Saunders  (1974)  discussed  several  other  estimators  of  3  .  First, 

,  U 

the  method  of  moments  estimator  6*  *  (rs)  was  shown  to  be  a  consistent 

w 

estimator  of  3  .  The  corresponding  mme  of  o  was  o*  ■  (s/B*  +  B*/r  -2) 

In  addition,  the  geometric  mean  $  *  (I  X  )^n  was  shown  to  a  consistent 

i-1  x 

estimator  of  B  as  well  as  the  median  estimator  defined  by 


X(k  +  1)  lf 
(X(k)X(k+l)),S 


n  -  2k  +  1 
if  n  -  2k  , 


where  X^  denotes  the  k**1  order  statistic  of  the  sample.  Thus, 
the  corresponding  estimators  of  (1.2),  R(t  ;  a*,  B*)  ,  etc.,  ere  con¬ 
sistent. 

It  is  known  (Blrnbeum  and  Saunders,  1969b  and  Saundera,  1974)  that 
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for  null  a  (a  <  1)  the  mle  0  and  the  mme  6*  virtually  agree. 
However,  as  a  result  of  the  simulation  in  Section  4  of  this  paper  (see 
Table  2),  for  a  i  1  the  mean  squared  error  of  0*  is  much  larger  than 
that  of  0  for  small  samples. 

3.  BAYES  ESTIMATION  OF  RELIABILITY  FOR  KNOWN  MEDIAN  LIFE 
Assuming  that  the  scale  parameter  and  median  0  is  known,  Bayes 
estimators  of  the  reliability  function  R(t  ;  o,0)  are  obtained  in 
this  section  with  respect  to  squared  error  loss.  First,  a  proper 
conjugate  family  of  priors  is  used,  and  then  a  vague  prior  of  Jeffreys 
is  considered.  A  Bayesian  analysis  with  both  parameters  a  and  0  un¬ 
known  seems  to  be  mathematically  intractable.  Thus,  a  modified  Bayes 
estimator  of  R(t  ;  u,0)  is  proposed  for  the  vague  priors,  where  0  is 
estimated  by  the  moment  estimators  in  Section  2,  and  is  compared  with  the 
mle  by  Monte  Carlo  simulations  in  Section  4. 

The  likelihood  of  the  random  aample  Jtj ,  ....  for  known  0  is 

L(a|xj,  ....  xn>  ■  Kj  o’n  exp(-K2/a2)  ,  (3.1) 

trtiere 

K,  -  (2/27  0)“n  x  xj  -  02 

1  ±m\  •  — -  -  ■  ■ 

xiC(xi./0),‘  -  (B/x^)**] 

and 

n 

K7  -  E  (x./0  +  0/x.  -  2)/2  . 

1  i-1  1  1 
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As  a  prior  distribution,  let  a  have  an  inverted  gamma  distribution 
with  parameters  (v,5)  ,  denoted  IG(v,6)  ,  and  pdf 

g(o2)  •  [6  T  (v)]  *(6/o^)V  +  *  exp(-6/a^)  ,  >  0  ,  (v,6  >  0)  , 

where  T(*)  denotes  the  usual  gamma  function.  That  is,  o  ^  is  the 
square-root  of  a  gamma  distributed  random  variable.  After  straightforward 

integration  of  (3.1)  multiplied  by  g(a^)  ,  the  posterior  distribution  of 

2  —1 
a  ,  given  the  data,  is  IG(v  +  n/2  ,  6  +  K2>  .  Therefore,  a  t 

given  Xj,  .  xn  ,  is  distributed  as  the  square-root  of  a  gamma  random 

variable  with  parameters  (v  +  n/2  ,  (6  +  •  Thus,  Leona  1  of 

Padgett  and  Wei  (1977)  may  be  applied  to  obtain  the  Bayes  estimate  of 

R(t  ;  a,B)  with  respect  to  a  squared  error  loss  function  as 

*nM  m  . x 

1  1  9  n 

-  PCT2v  +  n  <  -  c  ((v  +  n/2)/(«  +  K2)fy  ,  (3.2) 

where  is  a  random  variable  having  student's  t  distribution  with  k 

degrees  of  freedom. 

If  there  is  very  little  information  available  concerning  a  ,  a 
vague  or  noninformatlve  prior  may  be  used  to  obtain  a  (proper)  posterior 
distribution  for  o  .  Since  the  transformed  random  variable 

Is  U 

*  (X^/0)  -  (0/X^)  for  known  0  has  a  normal  distribution  with 

2 

me An  aero  and  variance  a  for  each  i-1,  ...,  n  ,  the  Fisher's  inf or- 

2 

nation  In(a)  can  found  to  be  -n/o  .  Thus,  Jeffrey's  noninformatlve 
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1973)*  It  should  be  noted  that  for  the  case  that  both  a  and  6  are 
unknown  Box  and  Tlao's  (1973)  vague  prior  idea  of  taking  the  prior  as 
p(B|o)«  constant  and  p(a)<*  a  *  apparently  does  not  result  in  a 
mathematically  tractable  posterior  distribution. 

For  the  improper  prior  p(ot)«  a  *  ,  the  (proper)  posterior  pdf  of 
a  ,  given  Xj,  ..m  x  »  is 

)  .  2K2n/2  e»p(-K2/a2)  >  „  „  „  > 

r (n/2)  an  +  1 


f(a|i 


n 


where  is  defined  as  before.  Therefore,  the  Bayes  estimator  of 

R(t  ;  a,B)  with  respect  to  a  squared-error  loss  function  is  given  by  the 

expectation 


Vfc) 


J”*(-c/a)  f(o|x1,  ....  xn)  d  a 

E0  [$(-c  0'S]  , 


(3.3) 


Ik  U 

where  c  -  (t/B)  -  (6/t)  and  0  has  a  gamma  distribution  with 
parameters  (n/2  ,  .  Applying  Lemma  1  of  Padgett  and  Wei  (1977) 
to  the  expectation  in  (3.3)  gives  the  Bayes  estimator  as 


R^t)  -  P[Tn  <  -  c  (n/2K2>!sl  ,  t  >  0  ,  (3.4) 

where  T  has  student's  t  distribution  with  n  degrees  of  freedom, 
n 

In  the  case  that  the  median  lifetime  B  is  also  unknown,  but  some 
prior  information  about  o  is  known,  a  "modified  Bayes"  estimator  may 
be  obtained  from  (3.2)  by  choosing  appropriate  inverted  gamma  prior 
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parameters  and  using  one  of  the  estimators  of  B  given  in  Section  2. 

Also,  in  the  case  of  unknown  B  with  vague  prior  information  about  a 
a  modified  Bayes  estimator  may  be  obtained  from  (3.4)  by  replacing  B 

A 

with  B  ,  3*  ,  or  another  estimator  from  Section  2.  In  the  next 

section,  the  behavior  of  this  second  modified  Bayes  estimator  is  indicated 
by  Monte  Carlo  simulation  results  when  8*  is  used  for  8  in  (3.4).  De¬ 
note  this  estimator  by  R*(t)  •  The  simulations  indicate  that  R*(t)  is 
generally  better  than  the  mme  R*(t)  *  R(t  ;  a*  ,  6*)  for  all  a  and 

A  A  A 

as  good  as  the  mle  R(t)  *  R(t  ;  a  ,  B)  only  when  a  s  1  .  However, 
the  drawback  of  the  mle  is  that  B  must  be  approximated  iteratively  as 
discussed  in  Section  2. 

4.  SMALL  SAMPLE  COMPARISONS  OF  ESTIMATORS  OF  RELIABILITY 

For  small  samples  direct  analytical  comparisons  of  the  mle  R(t)  , 

R*(t),  and  modified  Bayes  estimator  R*(t)  given  in  Sections  2  and  3  are  not 
feasible  due  to  the  mathematical  complexity  of  these  estimators.  Hence, 

Monte  Carlo  simulations  were  performed  to  compare  the  biases  and  mean  squared 
errors  of  the  estimators.  For  several  values  of  the  parameters  a  and  6  > 
1000  random  samples  of  size  n(*  10,  30)  were  generated  from  (1.1).  For 

A  ^  A 

each  sample  the  mle  o  ,  B  ,  and  R(t)  ,  the  mme  a*  ,  B*  ,  and 
R*(t)  •  and  the  modified  Bayes  estimate  R*(t)  were  computed  for  several 
values  of  time  t  .  The  average  squared  error  (ASE)  and  average  estimated 
reliability  (AKR)  for  each  of  the  three  estimates  of  R(t)  were 
computed  over  the  1000  samples.  Table  1  shows  some  of  the  simulation 
results.  Since  6  is  a  scale  parameter  the  results  in  Table  1  do  not 
depend  on  the  specific  value  of  S  .  Time  t  was  chosen  as  multiples 
of  B  .  The  same  results  were  obtained  when  2000  samples  were  used  to 
generate  the  ASE  and  AER  instead  of  1000. 


v?*j**r 
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TABLE  1.  Averages  and  ASE(xlO  )  of  Estimators  of  R(t;  a,  6) 


(times  t.  *  iB/2,  i  •  1,  2,  3,  4) 


(a)  n 


R(t;a,B  ) 


1.000 

0.500 


0.998 

0.500 

0.051 


0.921 

0.500 

0.207 


0.760 

0.500 

0.342 

0.240 


0.556 

0.500 

0.467 

0.444 


MLE  R(t) 


1.000  0.000 

0.488  2.011 


0.995  0.008 
0.488  2.001 
0.051  0.209 


0.924  0.382 

0.488  1.964 


MME  R*(t)  Mod.  Bayes  R^(t) 


1.000  0.000 

0.488  2.011 


0.995  0.008 
0.488  2.000 
0.051  0.209 


0.924  0.385 

0.488  1.964 


0.198  1.032  0.198  1.032 


0.776  1.140 
0.488  1.820 
0.335  1.437 
0.221  1.190 


0.566  0.542 
0.493  0.649 
0.466  0.505 


0.776  1.172 
0.489  1.841 
0.334  1.456 
0.220  1.212 


0.569  0.776 

0.493  0.869 


1.000  0.000 

0.489  1.893 


0.990  0.022 
0.489  1.883 
0.063  0.234 


0.911  0.391 
0.488  1.845 
0.208  0.949 


0.767  1.071 
0.489  1.733 
0.339  1.355 
0.229  1.104 


0.567  0.727 

0.493  0.823 
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TABLE  1.  (cont'd) 

(b)  n  -  30 


True 

MLE  R(t) 

MME  R* 

(t) 

Mod.  Bayes  R^(  t:) 

a 

Time 

_  .  j 

R(t;a ,3) 

AER 

ASE 

AER 

ASE 

AER 

ASE 

MM 

mm 

1 

0.002 

0.002 

m 

0.003 

0.25 

*2 

mSm 

0.533 

0.532 

m 

0.523 

C3 

0.051 

0.052 

0.081 

0.052 

0.081 

0.057 

0.086 

MB 

B 

0.136 

0.923 

0.136 

0.918 

0.136 

0.50 

C2 

1 

MM 

0.514 

0.501 

0.514 

0.501 

0.505 

"3 

0.207 

0.204 

0.348 

0.205 

0.350 

0.208 

0.339 

mm 

wmm 

0.765 

0.366 

mmm 

0.377 

mm 

0.364 

HSJb 

0.501 

0.444 

BE9 

0.455 

wEM 

0.447 

1.00 

MM 

mm 

0.339 

0.418 

0.339 

0.440 

0.340 

0.429 

K9 

m 

0.236 

0.376 

0.237 

0.392 

0.240 

0.380 

mm 

0.597 

0. 185 

0.597 

0.258 

0.596 

0.253 

IIb 

wMm 

0.501 

0. 160 

0.501 

0.226 

0.501 

0.222 

3.00 

1 

0.446 

0.444 

0.163 

0.444 

0.245 

0.444 

0.241 

H 

0.407 

0.403 

0.186 

0.404 

0.271 

0.405 

0.265 

MB 

mmm 

0.559 

0.078 

IB 

0.140 

0.558 

0.138 

5.00 

t2 

Bis 

0.501 

0.065 

mm 

0.119 

0.501 

0.117 

C3 

0.467 

0.466 

0.068 

0.466 

0.466 

0.130 

Table  2  shows  some  of  the  ASE's  for  the  mle  and  mme  of  a  and 
6  .  As  Indicated  in  Section  2,  the  estimator  B*  is  in  agreement 

A  A. 

with  g  for  ct  <  1  ,  but  differs  from  6  considerably  in  mean 

squared  error  for  a  £  1. 

As  the  simulation  results  of  Table  1  indicate,  the  modified  Bayes 
estimator  R*  (t)  is  better  in  mean  squared  error  than  the  moment 
estimator  R*(t)  for  all  values  of  a  •  It  also  is  no  worse  than  the 

mle  for  a  £  1  ,  although  the  performance  of  the  mle  is  better  than  that 
of  R*(t)  for  a  >  1  .  The  biases  of  all  three  estimators  are  similar 

for  all  values  of  a 

In  addition  to  using  the  moment  estimator  8*  to  modify  the  Bayes 
estimator  (3.4) ,  the  other  consistent  closed  form  estimators  of  8  » 

v  ^  % 1/n  ~ 

6  *  (  H  X.)  and  the  median  8  •  stated  in  Section  2  were  used, 

n  i„l  i  n 

These  produced  slightly  larger  ASE  than  R*(t)  did.  Also,  if  it  is 
known  that  a  is  near  a  specific  value,  one  of  the  proper  inverted 
gamma  priors  of  Section  3  could  be  used  to  produce  a  modified  Bayes 
estimator  from  (3.2).  This  would  generally  be  a  better  estimator  than 
R*(t)  . 

5.  EXAMPLE 

As  an  example,  some  of  the  fatigue  failure  data  of  Blrnbaum  and 

Saunders  (1958)  will  be  used.  Specifically,  the  failure  times  for  101 

strips  of  6061-T6  aluminum  sheeting  which  were  recorded  under  testing 

with  periodic  loading  at  18  cycles  per  second  and  maximum  stress  of 

31,000  psi  will  be  used  to  estimate  the  reliability  by  R(t)  ,  R*(t)  , 
and  R*(t)  .  This  data  is  known  to  follow  approximately  the  distribution 


Newton's  method  yielded  the  mle  of  6  as  6  ■  1336.37  ,  and  the 

mle  of  a  was  a  *  0.31032.  The  moment  estimator  of  6  was 

g*  *  1340.86  ,  and  a*  *  0.31034.  Table  3  gives  the  estimates  of 
reliability  R(t  ;  a,  8)  for  several  values  of  t  . 

TABLE  3.  RELIABILITY  ESTIMATION  FOR  ALUMINUM  SHEETING 


t 

600 

800 

1000 

1200 

1400 

1600 

1800 

2000 

R(t) 

0.9960 

0.9527 

0.8258 

0.6357 

0.4404 

0.2806 

0.1677 

0.0954 

R*(t) 

0.9961 

0.9538 

0.8286 

0.6398 

0.4447 

0.2843 

0.1705 

0.0973 

R*(t) 

0.9955 

0.9522 

0.8274 

0.6394 

0.4448 

0.2849 

0.1716 

0.0988 

6.  SUMMARY  AND  CONCLUSIONS 

For  the  case  that  the  median  life  (or  scale  parameter)  6  is  known, 
Bayes  estimators  of  the  reliability  function  for  the  Blrnbaum-Saunders 
model  have  been  presented  with  respect  to  squared  error  loss  functions  and 
a  proper  conjugate  family  of  prior s  (square  root  of  gamma)  or  a  vague 
prior.  If  both  parameters  a  and  $  are  unknown,  Bayes  solutions  for 
reliability  in  a  compact  form  seem  to  be  extremely  difficult.  Hence, 
an  appealing  modified  Bayes  estimator  R*(t)  is  proposed  in  which  the 
value  of  6  is  replaced  by  the  moment  estimate  $*  in  the  Bayes  solution 
(3.4).  This  modified  Bayes  estimator  is  easily  computed  from  Student's 
t  distribution  whereas  the  maximum  likelihood  estimator  of  reliability 

A 

must  be  computed  by  an  iteration  procedure  for  the  mle  B  of  8  The 


modified  Bayes  estimator  has  been  shown  by  computer  simulation  results  in 
Section  4  to  be  preferable  in  the  sense  of  mean  squared  error  to  the 
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method  of  moments  estimator  R*(t)  for  all  values  of  a  .  In  addition 
to  being  computationally  simpler  than  the  mle,  the  simulations  indicate 
that  R*(t)  is  as  good  as  the  mle  in  mean  squared  error  for  a  £  1  . 

The  biases  of  all  three  of  the  estimators  of  reliability  are  approximately 
the  same* 

It  should  also  be  remarked  that  minimum  variance  unbiased  estimators 
of  a  ,  3  ,  or  R(t  ;  a,  6)  are  not  available  due  to  the  difficulty  in 

obtaining  jointly  sufficient  complete  statistics  for  a  and  6  .  The 

pdf  (1.1)  is  not  of  the  exponential  form. 
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